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INTRODDCT ION 


A stochastic spatial computer model addressing coastal resource 
problems In Louisiana is being refined and validated using thematic mapper 
(TM) imagery. Investigators at the Southeast Fisheries Center, National 
Marine Fisheries Service (NMFS), in Miami and at Louisiana State 
University's Center for Wetland Resources and Center for Energy Studies in 
Baton Rouge are participating in the research, which consists of two major 
activities: 

1. Measuring land (or emergent vegetation) and water and the 
length of the interface between land and water in TM imagery of 
selected coastal wetlands (sample marshes) 

2. Comparing spatial patterns of land and water in the sample 
marshes of the imagery to those in marshes simulated by the com- 
puter model 

We completed a significant amount of work during the fourth semiannual 
period. This report discusses our accomplishments in modeling by the NMFS 
group and in digital processing of TM data by the LSU group. Planned 
future work is briefly discussed in a final section. 

LSU activities during the fourth semiannual period focused on prepro- 
cessing TM images of brackish marsh sites and tabulating data on spatial 
parameters from TM images of the salt marsh sites. We continued to use the 
Fisheries Image Processing System (FIPS) maintained by NMFS in Slidell, 
Louisiana, to analyze the TM scene. FIPS uses a Sperry-Univac V77/600 com- 
puter, color image display device, and other hardware to process remotely 
sensed digital imagery. The software is a modified version of the Earth 
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Resources Laboratory Applications Software (ELAS) (Graham et al. 1984). 

NMFS activltes were concentrated on improving the structure of the 
model and developing a structure and methodology for calibrating the model 
with spatial-pattern data from the TM imagery. 


PREPROCESSING OF TM IMAGERY 

The study sites were selected during the third semester of the study 
(see Browder et al. 1986). They are located in salt and brackish marsh 
areas on two abandoned delta lobes of the Mississippi River located in 
southern Louisiana. The early Lafourche lobe was an actively progradlng 
delta within the last 1,800 years; the late Lafourche lobe was active as a 
main distributary of the river within the last 600 years. 

We selected salt marsh study sites on each delta lobe that correspond 
to the boundaries of five contiguous U.S. Geological Survey 7.5-min 
topographic maps. We then divided the individual topographic maps into 
four contiguous quarters, each encompassing an area 192 elements wide and 
192 scan lines long on the TM image. The intersection of the four quarters 
was aligned to correspond to the center point of each topographic map. 

Thus, a total of 40 salt marsh sites, 20 on each delta lobe, were available 
for the analysis. 

TM images of the brackish marsh sites were georeferenced to fit a 
Universal Transverse Mercator projection with a north-south orientation. 

The registration technique was identical to the method used to rotate TM 
Images of the salt marsh sites; the ELAS modules PMGC and PMGE were used to 
accumulate ground-control points, generate polynomial least-squares mapping 
equations, and resample the image using the bilinear interpolation tech- 
nique. Registration accuracies for the early and late Lafourche sites were 
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35 and 22 m, respectively. 

We will have fewer brackish narsh sample sites than salt marsh sites 
because of cloud cover, the upper boundary of the TM scene, the pre- 
sence of upland areas, and the Irregular boundaries of brackish marsh 
areas. We have estimated that approximately 15-25 brackish marsh sites 
will be available for the analysis. 

TECHNIQUE FOR SEGMENTING TM IMAGES INTO LAND AND WATER 

We continued to explore Image enhancement and gray-level thresholding 
techniques for segmenting the TM images into land and water. We visually 
evalated land and water discrimination In experimental Images representing 
three enhancement techniques: (1) a band-5 Image with the shorelines 

enhanced by high-pass filtering (Molk 1980), (2) a product image generated 
from bands 4 and 5, and (3) a ratio image derived from bands 3 and 5. The 
experimental images included the Leevllle oil field, an oil and gas field 
comprising numerous small and linear water bodies such as ponds, lakes, and 
man-made canals. 

The product and ratio images seemed to provide slightly better 
discrimination of small water bodies than did the edge-enhanced band-5 image. 
The product and ratio images were approximately equal in their ability to 
discriminate shorelines In the test area; the product image, however, had 
an advantage over the ratio image because the effects of multiplying 
bands 4 and 5 seemed to reduce the number of noisy pixels in open-water 
areas and thus eliminated the need to digitally made these areas. A pro- 
duct image including all of the study sites was then generated and 
rescaled for a 0-255 display. 
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We used the entroplc thresholding technique described by Pun (1981) 
to segment the product Image Into land and water pixels. This automated 
technique takes Into account the asymmetrical shape of the gray-level 
histogram to define a threshold between land and water pixels. The 
threshold is derived from an anisotropic coefficient, which is related to 
the geometric shape of the gray-level histogram. Pun's method was selected 
Instead of other thresholding techniques (e.g., Kirby and Rosenfleld 1979; 
Weszka and Rosenfleld 1978) because It was simple to Implement on our 
system and seemed to provide a reasonable depiction of the land-water ratio 
in the TM Images of the study sites. 

TABULATION OF SPATIAL PARAMETER DATA FROM THE TM IMAGES 

A total of 40 binary land-water images were generated from the product 
image of the salt marsh sites. These images were labeled and stored on 
tape. Algorithms based on sequential ELAS commands were set up for batch 
processing to generate the data for each of these sites for the spatial 
model. We measured the following spatial-pattern parameters: (1) total 
numbers of land and water pixels; (2) total numbers of water pixels by scan 
line and element column; (3) total number of land-water joins; (4) total 
number of water pixel sides adjacent to one, two, three, and four other 
water pixels; and (5) number of pixels by water body size. A significant 
amount of processing was completed for the salt marsh study sites during 
this semester of activity (Table 1). 

Table 2 lists the percentage of each area that is open water and the 
number of land-water joins in each of the 40 salt marsh study sites. The 
percentage of an area that is open water indicates the level of disin- 
tegration (i.e., an area that is 40 percent water is disintegrated to the 
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40 percent level). The number of land-water joins is synonymous with 
Interface and, multiplied by 2, gives the cross-product statistic. A plot 
of land-water joins, or interface, is given in Figure 1. Two aspects of 
this plot are wrth noting. First, data are missing at the lower end of 
the scale of marsh disintegration. Fortunately, we expect to obtain data 
In’ this range from the brackish marsh sites, which have yet to be pro- 
cessed. Second, the plot suggests that the land-water interface con- 
tinuously declines with disintegration levels greater than 40 percent. We 
think that several parabolic curves of interface versus disintegration 
level may exist in this data, but marsh simulations with models calibrated 
to these data will be necessary to define the trajectories. 

TECHNIQUE TO GENERATE SPATIAL AUTOCORRELATION STATISTICS 

Testing was completed on some new ELAS software developed under 
contract by Delta Data Systems, Inc., of Picayune, Mississippi. The soft- 
ware computes the general cross-product statistic (R), the expected value 
and variance of R, and a test statistic based on a z-approxlmation from 
each binary land-water image derived from the TM data. We are using these 
parameters to quantify the level of spatial autocorrelation present in the 
distribution of water bodies in each study site. The software utilizes 
equations from Upton and Fingleton (1985). 

We encountered some Some operational problems with the new software. 
First, some elaborate modifications were required to accommodate the 16-bit 
architecture and limited memory of the host computer. These problems, 
coupled with the complexity of the task, resulted in lengthy execution 
times. An average of seven days was required to complete the computations 
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for a single 192-X-192 Image. Second, the FIPS computer Is not stable 
enough to complete batch jobs requiring 24 hours or more of execution time. 
We estimated that about one calendar year mould be required to complete the 
tabulations for the 40 salt marsh study sites, given the lengthy execution 
times, other demands placed on the system, and the Instability of the host 
computer. 

Fortunately, an alternative mas available. Upton and Fingleton (1985) 
provide a second set of equations for calculating R, E(R), and var(R) by 
means of the number of land pixels, the number of mater pixels, and the 
number of land-water joins (which is equivalent to the length of the land- 
water interface) of a rectangular or square matrix. These equations are 
inappropriate for areas of other shapes or those having irregular boun- 
daries; however, they apply perfectly to our square sites. By switching to 
new algorithms based on these equations, we were able to keep spatial- 
pattern computations within the range of practicality. 

MODEL REFINEMENT AND CALIBRATION 

The model has two weighting factors that affect the probability that a 

given pixel will convert to water at the next iteration and that, there- 
fore, affect the pattern of land and water In the simulated marsh. The 

first weighting fac tor, W, is related to the number of sides on which each 

land pixel is bor dered by mater pixels. The second probability factor, G, 
relates to the position of each pixel in the marsh relative to the main 
water body (in this case, the Gulf of Mexico). As the model originally 
was Vrltten, even when G > 0, only those water pixels Initially bordering 
the main mater body had heightened probabilities of converting to water. 

In testing the sen sltlvlty of this model to the two weighting factors, we 
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found that G did not appreciably affect the cross-product statistic (R), 
which we plan to use to calibrate the aodel to simulate specific marshes 
from TM imagery. Essentially, we had only one functional weighting fac- 
tor, W, which might have provided too little latitude for fitting the land- 
water patterns of our simulated marshes to those of actual marshes. 

The parameter R was «ich more sensitive to variation in G in a second 
version of the model completed during this semester of activity. In this 
model, if G > 0, pixels throughout the marsh have a higher probability of 
converting to water if they border on the main water body at any time 
during the simulation. In sensitivity tests of this model, we found a 
strong and highly nonlinear relationship of R to G and W. This new ver- 
sion of the model seems to be much more suitable to our needs; although, 
owing to the nonlinear effect of W and G on R, the tad; of setting W and G 
on the basis of R will be more difficult than anticipated. Figure 2 is not 
sufficiently detailed to provide a thorough representation, but roughly 
depicts how R varies with respect to G and W in a marsh of 48 x 48 pixels. 
Each point represents the mean from five simulation runs. Since the execu- 
tion of even one simulation is fairly slow (see our third semi-annual 
progress report, Browder et al. 1986), the process of developing such a 
graph takes several weeks. Because R is dependent upon the size of the 
marsh and we plan to be working with 192-X-192 samples from the TM imagery, 
we decided to start at once on a graph of this sort for the larger marsh 
rather than spending any more time on the graph for the smaller marsh. The 
graph being prepared will be more detailed and, formatted as a look-up 
table, will help us select the weighting factors to best approximate the 
land-water patterns of marshes in the TM imagery. 
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Using digitized Information from a habitat map (Wicker et al 1980), 
the LSU group tabulated data on spatial-pattern parameters for four 96-x 
96-pixel areas within the late Lafourche salt marsh study site. The habi- 
tat maps are stored as part of the Map Overlay and Statistical System 
(MOSS) data base maintained by the U.S. Fish and Wildlife Service. (MOSS 
digital data converted from a polygon to a cell format were available for 
our use.) The data were for four subquarters of the northwestern quarter 
of the Mink Bayou quadrangle. The parameters tabulated were number of land 
pixels, number of water pixels, number of land-water joins (J), and number 
of water pixels with sides adjacent to 1 (ADJ-1), 2 (ADJ-2), 3 (ADJ-3), and 
4 (ADJ-4) water pixels (Table 3). These were the first data on actual 
marshes that we had seen, and we used it to gain perspective on fine tuning 
and calibrating the model before we received any data from the TM imagery. 

The four marshes were in slightly different stages of disintegration, 
ranging from 27 to 57 percent water. For the other spatial-pattern parame- 
ters to be exactly comparable, the marshes would have to be at the same 
disintegration level. Keeping this in mind, we still thought the data 
suggested some independent variation of the side-adjacency data and J. 

For instance, the proportion of ADJ-4 pixels seems much higher in the 
southeastern quarter than in the northeastern quarter and the proportion of 
3-ADJ pixels seems much lower in the southeastern quarter than in the 
northeastern quarter, although the level of disintegration in these two 
marsh areas is about the same. This suggests that we might need more than 
one spatial-pattern parameter to set our weighting factors. Furthermore, 
It confirms our suspicion that two weighting factors should have separate 
effects. Clearly, the second version of the model is more suitable than 
the first one. 
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To further explore the relationship between J and ADJ-1, 2, 3, and 4, 
we executed 25 simulations (each with five replications) of a 48-X-48 pixel- 
marsh and, at the 50 percent level of disintegration, collected data on J 
and ADJ-all. Linear regressions (Table 4) indicated that J and ADJ-4 were 
highly correlated (R 2 - 0.89) and that over 40 percent of the variation in 
each of these variables could be explained by V, G, and their product. 

The numbers of water pixels with 2 and 3 sides adjacent to water were signi- 
ficantly related to only the G weighting factor and, therefore, appear 
Independent of the W weighting factor. The R 2 of the relationship between 
ADJ-3 and G was higher than that between ADJ-2 and G. All four side- 
adjacency variables were correlated with Interface (Table 5), but the rela- 
tionship between ADJ-3 and interface was the weakest (R 2 was only 0.27 for 
ADJ-3, whereas it was 0.95 for ADJ-1, 0.85 for ADJ-2, and 0.89 for ADJ-4. 

Several results of these analyses suggested that ADJ-3 would be a good 
variable to use in combination with R to select the best W and G to simu- 
late the spatial patterns of a given marsh. First, ADJ-3 was significantly 
related to G but not to W, and ADJ-3 had a stronger relationship with 
G than did ADJ-2. Secondly, ADJ-3 's correlation with R was weaker than 
that of the other side-adjacency variables; thus, W and G would be selected 
on the basis of two relatively independent variables. 

Tentatively, we plan to select starting values of W and G by deve- 
loping lsopleths of ADJ-3 and a corresponding look-up table to use in con- 
junction with the R-value look-up table to make our initial W and G selec- 
tions. By superimposing the two sets of lsopleths and determining the 
point. or general area of W and G where contours of the desired R and the 
desired ADJ-3 Intersect or most nearly intersect, we can select the W and G 
that will best approximate the spatial pattern of land and water in the 
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sample. It will be necessary to prepare pairs of Isopleths (one for R and 
one for ADJ-3) for five levels of disintegration (0.1, 0.2, 0.3, 0.4, and 
0.5) and to Interpolate between them to match TM marshes at all possible 
disintegration levels. Once we have found the general region of W and G by 
means of the Isopleths, we plan to approach the optimal values through an 
structured Iterative decision-making process. 


FUTURE ACTIVITIES 

During the upcoming semester we will complete all of the tabulations 
from TM Images of the salt marsh study sites and continue processing the 
brackish marsh study sites. We expect to complete the tabulations of spa- 
tial-pattern data on the brackish marsh sites by the end of the semester. 

Within the next few weeks we will have collected the necessary data 
from model simulations needed to construct the Isopleths and look-up tables 
for R and ADJ-3. The next step will be to develop the decision-making pro- 
cedure to fine-tune our estimates of weighting factors W and G. Then we 
will begin to select the W* s and G's for the salt marsh study sites. 
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1. Plot of the nunber of land-water joins versus the percentage of open- 
water area derived from the TM images of the AO salt marsh sample 
sites. 

2. Isopleth plot of the land-water interface as a function of G and W, 
constructed from mean values of five simulation runs for a 

4 8-x-4 8-pixel marsh. 
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Figure 1. Plot of the number of land-water joins versus the percentage of open-water area derived 
from the TM Images of the 40 salt marsh sample sites. 





Table 1* Progress to date In processing TM imagery (X's Indicate completed units) 
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Table 2. Percentage of open water area and number of land-water joins 
tabulated from TM images of the salt-marsh study sites. 


Number of 

Percentage land-water . 

Quadrangle name Qiarter water joins 


Leeville 


Mink Bayou 


Caminada Pass 


Bay Tambour 


Pelican Pass 


Grand Bayou du Large 


Lake La Graisse 


Central Isles Demieres 


Cocodrie 


NW 

42.939 

9,093 

NE 

34.926 

10,025 

SE 

39.912 

10,118 

SW 

46.742 

6,115 

NW 

26.590 

7,166 

NE 

24.536 

5,519 

SE 

29.435 

8,420 

SW 

29.574 

7,937 

NW 

62.093 

6,778 

NE 

82.867 

1,341 

SE 

99.265 

101 

SW 

58.019 

3,884 

NW 

51.180 

5,818 

NE 

76.546 

3,047 

SE 

88.723 

2,245 

SW 

84.378 

3,515 

NW 

98.964 

341 

NE 

82.598 

2,433 

SE 

86.683 

2,100 

SW 

97.244 

385 

NW 

53.833 

3,431 

NE 

67.255 

3,621 

SE 

40.435 

7,174 

SW 

83.716 

1,111 

NW 

91.132 

2,079 

NE 

99.278 

292 

SE 

100.000 

0 

SW 

99.997 

4 

NW 

75.770 

3,439 

NE 

90.053 

1,605 

SE 

95.182 

748 

SW 

93.533 

745 

NW 

33.434 

7,430 

NE 

68.894 

6,457 

SE 

85.786 

3,681 

SW 

54.715 

9,692 

NW 

62.972 

4,612 

NE 

30.735 

8,114 

SE 

42.377 

7,774 

SW 

36.396 

6,482 


Dog Lake 




Table 3. Spatial—pattem parameters tabulated from four subquarters of the northwestern quarter of 
the 1978 Mink Bayou habitat map (Wicker et al. 1980). 



.Table 4. Results of the multiple regression of selected spatial-pattern parameters 
from 25 simulations of a 48-x-48-pixel marsh on probability weighting 
factors of the model. 



Regression coefficients of 
the independent variables 



• . 


Dependent 

variable 

W 

G 

W X G 

Constant 

r2 8 

SE of 
estimate 

Sig. F 

L-W Joins 
SE 

Sig. T 

-1.09141 

0.44948 

(0.0242) 

-1.53730 

0.44948 

(0.0026) 

0.00313 

0.00172 

(0.0826) 

983.59248 

117.55361 

(0.0000) 

0.41422 

0.33053 

309.35 

0.0094 

2-ADJ 

SE 

Slg. T 


-0.16651 

0.06933 

(0.0248) 


137.7444 

18.13118 

(0.0000) 

0.2005 

0.16577 

65.77 

0.0248 

3-ADJ 

SE 

Sig. T 


-0.24769 

0.06531 

(0.0009) 


217.98880 

17.07962 

(0.0000) 

0.38479 

0.35804 

61.95 

0.0009 

4-ADJ 

SE 

Sig. T 

0.41324 

0.20091 

(0.0523) 

0.76386 

0.20091 

(0.0010) 

-0.00149 

0.0007682 

(0.0655) 

559.74664 

52.54419 

(0.0000) 

0.43700 

0.35657 

138.27 

0.0063 


8 The second value In each set of r 2 is adjusted R^. 


Note: 


L-w joins - number of land-water joins. 
ADJ-1 ■ number of water pixels bordered on 

ADJ-2 * number of water pixels bordered on 

ADJ-3 * number of water pixels bordered on 

ADJ-4 - number of water pixels bordered on 


one side by another water pixels, 
two sides by other water pixels, 
three sides by other water pixels, 
four sides by other water pixels. 


.Table 5. Results of a simple linear regression of the side-adjacency variables on the 
land-water Interface. 


Dependent 

variable 

Regression coefficient 
(interface) 

Constant 

r2 8 

SE of 
estimate 

Sig. F 

ADJ-1 

0.12221 

-25.97183 

0.95351 

10.42236 

0.0000 

SE 

0.00563 

4.02513 

(0.95149) 



Sig. T 

(0.0000) 

(0.0000) 




ADJ-2 

0.17564 

0.29345 

0.85045 

28.44528 

0.0000 

SE 

0.01536 

10.98562 

(0.84395) 



Sig. T 

(0.0000) 





ADJ-3 

0.10706 

107.89119 

0.27401 

67.30131 

0.0072 

SE 

0.03634 

25.99189 

(0.24245) 



Sig. T 

(0.0072) 

(0.0004) 




ADJ-4 

-0.43045 

986.66296 

0.89131 

58.05365 

0.0000 

SE 

0.03134 

22.42042 

(0.88658) 



Sig. T 

(0.0000) 

(0.0000) 





8 The second value in each set of is adjusted R^. 

Note: L-W joins * number of land-water joins. 

ADJ-1 “ number of water pixels bordered on one side by another water pixel. 

ADJ-2 - number of water pixels bordered on two sides by other water pixels. 

ADJ-3 * number of water pixels bordered on three sides by other water pixels. 

ADJ-4 * number of water pixels bordered on four sides by other water pixels. 




